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Abstract 

In our recent work, we proposed the design of perfect reconstruction orthogonal wavelet filterbanks, called gmph- 
QMF, for arbitrary undirected weighted graphs. In that formulation we first designed "one-dimensional" two-channel 
filterbanks on bipartite graphs, and then extended them to "multi-dimensional" separable two-channel filterbanks for 
arbitrary graphs via a bipartite subgraph decomposition. We specifically designed wavelet filters based on the spectral 
decomposition of the graph, and stated necessary and sufficient conditions for a two-channel graph filter-bank on 
bipartite graphs to provide aliasing-cancellation, perfect reconstruction and orthogonal set of basis (orthogonality). 
While, the exact graph-QMF designs satisfy all the above conditions, they are not exactly fc-hop localized on the 
graph. In this paper, we relax the condition of orthogonality to design a biorthogonal pair of graph-wavelets that can 
have compact spatial spread and still satisfy the perfect reconstruction conditions. The design is analogous to the 
standard Cohen-Daubechies-Feauveau's (CDF) construction of factorizing a maximally-flat Daubechies half-band 
filter. Preliminary results demonstrate that the proposed filterbanks can be useful for both standard signal processing 
applications as well as for signals defined on arbitrary graphs. 



Note: Code examples from this paper are available at http://biron.usc.edu/wiki/index.php/Graph_FilterbankSi 



EDICS Category: DSP-WAVL, DSP-BANK, DSP-MULT, DSP-APPL, MLT 

I. Introduction 

A. Motivation 

Graphs provide a flexible model for representing data in many domains such as networks, computer vision, and 
high dimensional data-clouds. The data on these graphs can be visualized as a finite collection of samples, leading 
to a graph-signal, which can be defined as the information attached to each node (scalar or vector values mapped 
to the set of vertices/links) of the graph. Examples include measured values by sensor network nodes or traffic 
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measurements on the links of an Internet graph. The formulation of datasets as graph-signals has been subject 
to a lot of study recently, especially for the purpose of analysis, compression and storage. Major challenges are 
posed by the size of these datasets, making them difficult to visualize, process or analyze. This has led to a recent 
interest in extending wavelet techniques to signals defined on graphs. These techniques can provide multiresolution 
representations, so that smaller graphs with smooth approximations of the original graph-signals can be obtained 
and used for processing. Moreover, these signal representations can be used to develop local analysis tools, so that 
a graph-signal can be processed "locally" around a node or vertex, using data from a small neighborhood of nodes 
around the given node. 

The trade-off between spatial and frequency localization is fundamental in the study of wavelet transforms for 
regular signals. This trade-off is being studied for graph signals as well HI, Q, Q. We can say that a graph 
signal, for example an elementary basis function in a wavelet representation, is localized in the vertex domain if 
most of the signal energy is concentrated in a k-hop localized neighborhood around a vertex, so that the degree 
of localization will depend on how large k is. Similarly, a graph signal can be said to be localized in the spectral 
domain, defined in terms of the eigenvalues and eigenvectors of the graph Laplacian matrix, if the projection of 
the signal onto these eigenvectors has most of its energy concentrated in a band of graph-frequencies around a 
center frequency. 

There has been a significant recent interest in the extension of wavelet transforms to graph signals, including 
wavelets on unweighted graphs for analyzing computer network traffic JU, diffusion wavelets and diffusion wavelet 
packets Q, lH, 0, the "top-down" wavelet construction of jSl, graph dependent basis functions for sensor network 
graphs f9l, lifting based wavelets on graphs ifTOl . ifTTl . ifTlll . multiscale wavelets on balanced trees lfT3l . spectral 
graph wavelets lUl, and our recent work on graph wavelet filterbanks ||2]|. 

In this paper, we focus on the problem of designing wavelet transforms that are invertible, compactly supported 
on the graph and critically sampled (CS). Critical sampling, i.e., the number of wavelet coefficients generated by the 
transform is equal to the number of vertices in the graph, is important in order to achieve a compact representation 
(e.g., for compression). Also, in a critically sampled transform a subset of vertices will include low frequency 
information. This leads to a natural approximation of the original graph by a smaller graph (containing only those 
vertices and their corresponding graph signal). 

Of the above mentioned approaches, only the lifting based approaches ifTOll . ifTTll . the wavelets on balanced trees 
|[T3l and the graph wavelet filterbanks ||2J achieve critical sampling, although in all cases there are some restrictions 
on the kinds of graphs on which the transforms can be defined. 

In terms of localization, the lifting-based approaches ifTOl . lITTI are compactly supported in the graph, i.e., their 
basis functions are strictly fc-hop localized the vertex domain, which implies that the output at each vertex n can 
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be computed exactly from the data at the vertex n and a A;-hop neighborhood around it. Instead, in our filterbank- 
based approach ||2l basis functions had good vertex domain localization but were not compactly supported. Compact 
support could be achieved by approximating their spectral response by a polynomial, but this comes at the expense 
of introducing a small reconstruction error (i.e., there was no longer perfect reconstruction). 

The main contribution in this paper is to extend the filterbank approach to design graph wavelets of 1*21 so as 
to achieve compact support. Note that the lifting based techniques lITOl . ifTTll are also compactly supported but 
they are vertex domain designs for which it is difficult to control the quality (e.g., localization) of their spectral 
representation. Instead, in the work presented here filters are designed in the spectral domain while guaranteeing 
compact vertex domain support, so that we can directly control explicitly the trade-off between localization in the 
vertex domain and the spectral domain. 

As in 111, the building blocks our design are two channel wavelet filterbanks on bipartite graphs, which provide 
a decomposition of any graph-signal into a lowpass (smooth) graph-signal, and a highpass (detail) graph-signal. For 
arbitrary graphs the filterbanks can be extended in two ways: a) by implementing the proposed wavelet filterbanks 
on a bipartite graph approximation of the original graph, which provides a "one-dimensional" analysis, or b) 
by decomposing the graph into multiple link-disjoint bipartite subgraphs, and applying the proposed filterbanks 
iteratively on each the subgraphs (or on some of them), leading to a "multi-dimensional" analysis Q. 

In we showed that downsampling/upsampling operations in bipartite graphs lead to a spectral folding 
phenomenon, which is analogous to aliasing in regular signal domain. We utilized this property to propose two 
channel critically sampled wavelet filterbanks, called graph-QMF, on arbitrary undirected weighted graphs. We 
specifically designed wavelet filters based on the spectral decomposition of the graph, and stated necessary and 
sufficient conditions for a two-channel graph filterbank on bipartite graphs to provide aliasing-cancellation, perfect 
reconstruction and orthogonal set of basis (orthogonality). While the exact graph-QMF designs satisfy all the above 
conditions, they are not compactly supported on the graph. In order to design compactly supported graph-QMF 
transforms, we performed a Chebychev polynomial approximation of the exact filters in the spectral domain, and 
this incurred error in the reconstruction of the signal, and loss of orthogonality. Here, we propose an alternative 
to graph-QMF design where we relax the conditions of orthogonality, and design a biorthogonal pair of graph- 
wavelets, which we call graphBior. These new designs lead to a representation that is exactly A;-hop localized, 
and still satisfies the perfect reconstruction conditions. This design is analogous to the standard Cohen-Daubechies- 
Feauveau's (CDF) [14] construction to obtain maximally half-band filters. Even though these filterbanks are not 
orthogonal, we show that they can be designed to nearly preserve energy. In particular, we compute expressions 
for Riesz bounds of the filterbanks, and choose graph-wavelets with the maximum ratio of lower and upper Riesz 
bounds. 
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We will show that it is possible to design these filterbanks based on both the normalized Laplacian and the 
random-walk Laplacian, leading to nonzeroDC graphBior and zeroDC graphBior designs, respectively. As the 
name suggests, the zeroDC design has the advantage of leading to highpass operators with a zero response for the 
all constant signal, which will be useful in applications in Euclidean space (e.g., when applying graph wavelets 
to regular domain signals such as images). Instead, in the nonzeroDC design the DC frequency corresponds to a 
signal where the value at each node depends on its degree. Finally, because our designs are biorthogonal, the norms 
of highpass and lowpass are not equal. For applications where a normalization is required, we propose unity gain 
compensation (GC) techniques for both types of designs. A comparison of our proposed design vis-a-vis existing 
transforms is shown in Table HI 



Method 


DC 


CS 


PR 


Comp 


OE 


GS 


Wang & Ramchandran |9| 


No 


No 


Yes 


Yes 


No 


No 


Crovella & Kolaczyk [4] 


Yes 


No 


No 


Yes 


No 


No 


Lifting Scheme IfTOl 


Yes 


Yes 


Yes 


Yes 


No 


Yes 


Wavelets on balanced trees 1131 


Yes 


Yes 


Yes 


Yes 


No 


Yes 


Diffusion Wavelets |5| 


Yes 


No 


Yes 


Yes 


Yes 


No 


Spectral Wavelets |JJ 


Yes 


No 


Yes 


Yes 


No 


No 


graph-QMF filterbanks (exact) |2| 


No 


Yes 


Yes 


No 


Yes 


No 


graph-QMF filterbanks (approx.) |2| 


No 


Yes 


No^ 


Yes 


No^ 


No 


nonzeroDC graphBior filterbanks 


No 


Yes 


Yes 


Yes 


No 


No 


zeroDC graphBior filterbanks 


Yes 


Yes 


Yes 


Yes 


No 


No 



TABLE \: Comparison of graph wavelet designs in terms of key properties: zero highpass response for constant 
graph-signal (DC), critical sampling (CS), perfect reconstruction (PR), compact support (Comp), orthogonal 

expansion (OE), requires graph simpUfication (GS). 



The outline of the rest of the paper is as follows: in Section|lT} we introduce some notations, a general formulation 
of two channel wavelet filterbank on graphs and the graph-QMF filterbanks proposed in [2J, which are orthogonal 
and perfect reconstruction. In Section |IIIJ we describe proposed nonzeroDC graphBior filterbanks on graphs, and 
in Section |IV] we design and describe the properties of zeroDC graphBior filterbanks. In Section [V| we describe 
extension of proposed filterbanks to arbitrary graphs via bipartite subgraph decomposition, and multiresolution 



implementation. In Section VI we conduct some experiments to demonstrate the properties and applications of the 



proposed filterbanks. Section VII concludes our paper 



II. Preliminaries 

A graph can be denoted as G = {V,E) with vertices (or nodes) in set V and links as tuples in E. The 

graphs considered in this work are undirected graphs without self-loops and without multiple links between nodes. 

The links can only have positive weights. The size of the graph N = \V\ is the number of nodes and the geodesic 

'The exact Graph-QMF solutions are perfect reconstruction and orthogonal, but they are not compact support. Localization is achieved 
with a matrix polynomial approximation of the original filters, which incur some loss of orthogonality and reconstruction error, which can 
be arbitrarily reduced by increasing the degree of approximation. 
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distance metric is given as dg{i,j), which represents the sum of link weights along the shortest path between 
nodes i and j, and is considered infinite if i and j are disconnected. Define Mj^n as the j-hop neighborhood of 
node n, i.e., Mj^n = {"i G V,dg(i,j) < j}. Denote the identity matrix as I, and let 5„ be an impulse function, 
i.e., (5„(n) = 1 and 6n{rn) = for all m ^ n. Define < f, g >= f^g as the inner-product of vector f and g, 
where (.)^ is the transpose operator. Define A = [wij], the adjacency matrix of the graph, di the degree (sum of 
link-weights) of node i, and D = diag{di]i=i^2,...N, the diagonal degree matrix of graph, so that L = D — A is 
the unnormalized Laplacian matrix of the graph. 

A. Graph vertex domain 

We define a signal / : V — > M on a graph as a set of scalars, where each scalar is assigned to one of the vertices 
of the graph. This can be extended to vector values at each vertex. Further, a graph based transform is defined as 
a linear transform T : — applied in the vertex domain, such that the operation at each node n is a linear 
combination of the value of the graph-signal f{n) at the node n and the values f{m) on nearby nodes m G Mj^n, 
i.e., 

y(n) =< T(n,.), f >= r(n, n)/(n) + T{n,m)f{m), (1) 

where T(n, .) is the n*'^ row of the transform T. In analogy to the 1-D regular case, we would sometimes refer to 
graph-transforms as graph-filters, and T(n, .) for n = 1, 2, ...N as the impulse response of the transform T at the 
n*^ node. Note that due to the irregularity of links, the impulse response of a graph filter varies from one vertex to 
the other. A desirable feature of graph filters is spatial localization, which typically means that the energy of each 
impulse response (i.e., each row) of the graph filter is concentrated in a local region around a node. In this paper, 
we use the definition proposed in ||3l to define the spatial spread of any signal f on a graph G as: 

A2,(f) :=inf ^ (2) 

M 112 

Here, |f | P}j=i,2,...,Af can be interpreted as a probability mass function (pmf) of signal f, and A|,(f) is 

the variance of the geodesic distance function dg{i,.) : V — )■ M at node i, in terms of this spatial pmf. Thus, 
A^(T(n, .)) should be small for all n = 1,2,. ..N for good spatial localization. In our analysis, we compute the 
spatial spread of a graph transform T to be the average of the spatial spread Q of impulse responses over all 
vertices, i.e., 

1 ^ 

^G(T):=^5;A2,(T(n,.)) (3) 

n=l 
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B. Graph spectral domain 

As in our previous design O, we use the symmetric normaUzed Laplacian matrix C = D~^/^LD~^/^ to 
define spectral properties of the graph. In this paper, we use the same Laplacian matrix C to design nonzeroDC 
filterbanks. Because C is a real symmetric matrix, it has a complete set of orthonormal eigenvectors, which we 
denote by {ui};=o,i,2,...,Af-i- These eigenvectors have associated real, non-negative eigenvalues {Xi}i=o,i,2,...,N-i 
satisfying Cui = A;U;, for / = 0, 1, 2, A^ — 1, and ordered as: Aq < Ai < ...Aat-i. In our analysis, we assume the 
graph to be connected Zero appears as a unique minimum eigenvalue of the graph and the maximum eigenvalue 
is always less than or equal to 2, with equality if and only if the graph is a bipartite graph I115i Section 2]. 

We denote the spectrum of the graph by (j{C) := {Aq, Ai, ...Aat-i}. A graph signal f is represented in the 
spectral domain as its projection onto the eigenvectors, denoted as : /(A/) =< f, >. An eigenspace Vx is 
defined as the space spanned by eigenvectors of C associated with eigenvalue A, and the eigenspace projection 
matrix is defined as: 



Pa := 



T 
I ) 



A,=A 

where is the transpose of eigenvector u;. The eigenspace projection matrices are idempotent and Pa and P-y 
are orthogonal if A and 7 are distinct eigenvalues of the Laplacian matrix, i.e., 

PaP^ = <5(A-7)Pa, (4) 

where 6{X) is the Kronecker delta function. Further the sum of all eigenspace projection matrices for any graphs is 
an identity matrix. Closely related to the symmetric Laplacian matrix C is the random-walk graph Laplacian, which 
is defined as Cr := D^^L. We use Cf to design zeroDC filterbanks. Note that Cr has the same set of eigenvalues 
as C, and if u/ is an eigenvector of C associated with A;, then D~^u; is an eigenvector of Cr associated with the 
eigenvalue A^. Similar to (|2]), the spectral spread of a graph signal c an be defined as: 



(5) 

I III Mo ' ' L J I 

where {[/(A)]2/||f||2}A=Ao,Ai,...,AN_i is the pmf of f across the spectrum of the Laplacian matrix, and jj, is the 
mean of A with respect to this pmf. A^(f) computes the variance of A with respect to the spectral pmf function. 
Thus, for a signal to have good spectral localization, the value of A^(T(n, .)) should be small for all vertices. In 
our analysis, we compute the spectral response of any transform T as the average spectral response of impulse 
^For graphs with multiple separate connected components, we analyze each component as a separate graph. 
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responses at all vertices, i.e., 

1 ^ 

mX)\':=j^^\T{n,.m', (6) 

i=l 

and use |T(A)p as the spectral pmf function to compute the spectral spread of ij^ 
C. Spectral graph filters 

For designing compact support wavelet filters on graphs we use the same approach as in [2], and define analysis 
wavelet filters Hj and Gj in terms of spectral kernels hi{X) and gi{X) for i = 0,1 respectively. The corresponding 
transform matrices are represented as: 



Xea{G) ^rj-^ 

These filters have the following interpretation: the output of a spectral filter with kernel h{X) can be expanded 
as in = Hf = ^^xerriG) ^i{X) Px^, where fx = P^f is the component of input signal f in the A-eigenspace. 
Thus, filter H either attenuates or enhances different harmonic components of the input signal depending upon the 
magnitude of h{X). Therefore, we will also refer to h(X) as the spectral response of filter H. For a general kernel 
function, the filtering operations corresponding to Hj and Gj may not have compact support, and would require 
a full spectral decomposition of Laplacian matrix. However it has been shown in [1], that the spectral response 
can be approximated as a polynomial of degree K, and the corresponding filters can be computed iteratively with 
K one-hop operations at each node. Further, any graph filters with a K degree polynomial spectral response are 
exactly K-hop localized (compact support) O] Lemma 5.2], and can be efficiently computed without diagonalizing 
the Laplacian matrix. The computational complexity of the filtering operations in the polynomial case, reduces to 
0{K\E\) for degree K and \E\ number of links in the graph. Thus, the degree K in case of polynomial spectral 
response can be interpreted as the length of the corresponding spectral filters |^ 

D. Spectral wavelet filterbanks 

In Q, we described the construction of a two-channel wavelet filterbank on a bipartite graph B = {L, H, Ef\ 
characterized by filtering operations {Hj, Gj}j=o,i and a function which provides a decomposition of graph- 

' Note that the definitions of spread presented here are heuristic and do not have a well-understood theoretical background. Another 
definition of spectral spread in graphs is given in f3l. If the graph is not regular, the choice of which Laplacian matrix {C or jC) to use for 
computing spectral spreads also affects the results. The purpose of these definitions and the subsequent examples is to show that a trade-off 
exists between spatial and spectral localization in graph wavelets. 

*Note that, having a polynomial spectral response for compact support is necessary only in case of spectral graph filters. There can be 
non-spectral graph-filters, for example, graph wavelets proposed in (4), that have compact support without being a polynomial in the spectral 
domain. 

bipartite graph G — (L, H,E) is a graph whose vertices can be divided into two disjoint sets L and H, such that every link connects 
a vertex in L to one in H. 
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signal f into a lowpass (approximation) graph-signal f;o«, and a highpass (details) graph-signal component fhigh- The 



analysis side 




synthesis side 




Gi 



^ high 



Fig. 1 : Block diagram of a two-channel wavelet filterbank; on graph. 

transforms Hj, for i = 0, 1 of the two channels are graph transforms with spectral kernels hi and respectively, 
as given in (|7]l. In the analysis side of the filterbank, the input signal is first operated upon by transform Hq in 
the lowpass channel and Hi in the high pass channel. A subsequent downsampling upsampling (DU) operation, 
discards and replaces with zeros the output coefficients on the set H in the lowpass channel and on the set L in 
the highpass channel. Since L and H are disjoint and complementary subsets of vertex set V, the retained set of 
output coefficients is critically sampled. Algebraically, the DU operation can be represented with a function /3(n), 
such that P{n) = 1, if node n £ L and /3(n) = —1 if node n £ H. Thus, the DU operation in the lowpass channel 
is given as |(1 + /3(n)) and in the highpass channel ^(1 — /3(n)). It can also be written in the matrix form as 
^(I + J^) for lowpass channel and ^(I — J^) for highpass channel, where = diag{P} is a diagonal matrix. 
The output signal i^u of the DU operation in the two channels is expressed in Figure |2] 



f = 



f = 




A 



A 



i 




in 



|(I + J/3)f 



= 1(1- J/3)f 



Fig. 2: DU operations in the proposed two-channel filterbank. 



The overall transfer matrix of the two channel filterbank can be written as: 



T = ^Go(I + J;3)Ho + ^Gi(I- J;3)Hi 



-(GqHo + GiHi) + -(G0J/3H0 — G1J/3H1), 



(8) 



where Teg is the transfer function of the filterbank without the DU operation and TaHas arises primarily due to 
the DU operations. Using the above formulation, we derived the following result^ 
'see 12 for proofs and details. 
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Theorem 2.1 ("Spectral folding phenomenon |T, Prop. 1]): Given a bipartite graph B = {L,H,E) with Lapla- 
cian matrix C and with the binary function (3 defined as above, if ux is the eigenvector of a unique eigenvalue A 
of graph B then 

J/3UA = ±U2_A. (9) 

The ambiguity of it sign appears in Q since both ux and —ux can be eigenvectors of B for eigenvalue A. Further, 
if the eigenvalue A appears with multiplicity greater than 1, and if Pa is the projection matrix corresponding to 
the eigenspace Vx, then 

JpPx = P2-AJ/3- (10) 



According to Q and ( [TO) ) the eigenvector (or eigenspace) of eigenvalue A in a bipartite graph changes to the 
eigenvector (or eigenspace) of eigenvalue 2 — A after multiplying with J^. In other words, the eigenspace folds 
across the imaginary axis at A = 1, hence we call this a spectral folding phenomenot^ 

Theorem 2.2 (Perfect reconstruction property |2, Sec. III.B]): Given a bipartite graph B = {L,H,E), the 
filtering operations {Hj,Gi} and the binary function /3(n) as defined above, a necessary and sufficient condition 
for the perfect reconstruction in the two channel filterbanks is that for all A in (j{B), 

go{X)ho{X) + giiX)hi{X) = 2, 
ffo(A)Ao(2 - A) - giiX)hi{2 - A) = 0. (11) 



The advantage of representing PR conditions as in (111 is that the filterbank can be designed in the spectral 



domain of the graph by designing spectral kernels which satisfy ( [TT] ). These kernels can be designed as continuous 



functions of A G [0 2], which obviates the need to evaluate (111 only at the spectrum cr{B), thus making the design 



independent of the structure of the graph. The graph-QMF solution provided in 121, satisfies (111, and is also 



orthogonal. While there exist many exact graph-QMF solutions, we show in Section III-A (see last paragraph) that 
none of the solutions has compact support. Therefore, we shift our focus to biorthogonal solutions which are PR 
and have compact support. 

III. NONZERODC GRAPHBIOR FILTERBANKS 
Given the parameter k, in order to design /c-hop localized filters to satisfy the perfect reconstruction conditions 



given in (111, we need to design four polynomial spectral kernels of degree k, namely lowpass analysis kernel 



^Note that spectral folding phenomenon only occurs if the binary function /3(n) is defined on one of the natural partitions L or _ff of the 
bipartite graph B — {L, H, E), and not for any other partition. 
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/io(A), highpass analysis kernel /ii(A), lowpass synthesis kernel go{X), and highpass synthesis kernel 51(A). If we 
choose analysis and synthesis highpass kernels to be: 

hi{X) = 50(2 -A) 

51(A) = /io(2-A), (12) 



then, ( [TT) reduces to a single constraint for all eigenvalues, given as: 

^o(A)5o(A) + ho{2 - \)go{2 - A) = 2. (13) 



Further, define ^(A) = ho{X)go{X), then (13 1 can be written as: 

p{X) + p{2 - X) = 2. (14) 



In our approach, we first design /io(A) and 50(A), and then hi{X) and 51(A) can be obtained using (12l. Further, 
since p{X) is the product of two lowpass kernels, it is also a lowpass kernel. Therefore, the objective is to design 
p{X) as a polynomial half-band keme^ which satisfies (13 1, and then obtain kernels ^o(A) and ^i(A) via spectral 



factorization. This design is similar to the design proposed by Cohen-Daubechies-Feauveau |14| for finding a 
maximally fiat pair of lowpass and highpass filters under the given length constraint, and then dividing up the 
residual factors between the two filters in a way that makes the basis function nearly orthogonal. The following 
result is useful in our analysis: 



Proposition I: If /io(A) and 50(A) are polynomial kernels, then any j5(A) = /io(A)5o(A), which satisfies (14i for 
all A G [0 2], is an odd degree polynomial. 



Proof: By changing the variable so that A = 1 + /, we can write ( 14 1 as 



p{l + l)+pil-l) = 2, (15) 

where p{l + /) = /io(l + 05o(l + /) and / G [—1 1]. If ho{l) and 5o(Z) are polynomial in I then the functions 
p{l + /) and p{l — I) are also polynomials in and can be expressed as: 

K 

p{l + l) = ^Ck{l)^ 
k=0 
K 

p{l-l) = ^Cki-l)". (16) 

k=0 

^An ideal half band kernel h{X) in the spectral domain of a bipartite graph can be defined as a kernel with h{\) = 1 for A < 1, and 
h{X) = otherwise. 
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Using ([T6]) in ([T5]), we get: 



K 



K/2 



(17) 



fc=0 



A:=l 



Thus p{l + ^) + p(l — I) is an even polynomial function of I. In order for both (15 1 and (17i to be true for all 



/ G [— 1 1], Co = 1 and all other even power coefficients c„ in the polynomial expansion of j5(l + /) must be 0. 
Therefore, the solution p(l + I), expressed as: 

K 



p{l + l) = l + ^C2n+l/ 



2n+l 



(18) 



n=0 



is an odd degree polynomial. Thus, ignoring the trivial case p(l + Z) = 1, the highest degree of p(l + 1) and p{l — l) 
(and hence p{\)) is always odd. ■ 

A. Designing half-band kernel p{\) 



The following known results help us prove the existence of a polynomial p{\) that satisfies (15 1, and obtain its 
spectral factorization: 

Lemma 1 (Bezout's identity / f76] prop. 3.13]): Given any two polynomials a{l) and b{l) of continuous variable 

I, 

a{l)x{l) + b{l)y{l) = c{l), (19) 

has a solution [x{l), y{l)\, if and only if gcd{a{l),h{l)) divides c(Z), where gcd{a{l),b{l)) refers to the greatest 
common divisor of polynomials a{l) and b{l). 

Theorem 3.1 (Complementary Filters / [76] prop. 3.13]): Given a polynomial kernel ho{l), there exists a comple- 



mentary polynomial kernel 5o(0 which satisfies the perfect reconstruction relation in ( 13 1, if and only if /io(l + /) 
and ho{l — I) are coprime. 



Proof: Let us denote a{l) = /io(l + I), b{l) = ho{l — I) = a{—l), and c(/) = 2. Then, (13 1 can be written in 

a{l)x{l) + a{-l)y{l) = 2, (20) 



the same form as (19 1, i.e.. 



We first note that if a polynomial solution [x{l), y{l)] of (20 1 exists, then 



a{-l)x{-l) + a{l)y{-l) = 2, 



(21) 



also has a polynomial solution. Subsequently, combining (20 1 and (21 1 and choosing ^o(l + = l/2(a^(0 + 
we find that 

a(-05o(l + + a(05o(l-0 = 2, (22) 
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also has a polynomial solution. However, based on Lemma[T| ( pO| ) has a polynomial solution if and only if gcd{ho{l+ 
/), /io(l — 0) divides c{l) = 2, which is a prime number. This implies gcd{hQ{l + /), ho{l — I)) is either 1 or 2 for 
all / G [—1 1], which is true iff /io(l + /) and ho{l — I) do not have any common roots. This implies that /io(l + 1) 
and ho{l — I) are coprime. ■ 
Corollary 3.2 ( I\j6\ exercise. 3.12]): There is always a complementary filter for the polynomial kernel (1 + Z)*^, 



I.e., 



{l + l)''R{l) + {l-l)''Ri-l) = 2 



(23) 



always has a real polynomial solution R{1) for k > 0. 

Proof: Let us denote a{l) = (1 + 6(0 = (1 - I)'', x{l) = R{1), y{l) = R{-1) and c{l) = 2. Then, 



can be written in the same form as (19 1. Since a{l) and b{l), in this case are coprime, therefore gcd{a{l), b{l)) = 1 



divides c{l) = 2. Hence, a polynomial R{1), which satisfies p3| ) always exists. ■ 
For a perfect reconstruction biorthogonal filterbank, we need to design a polynomial half-band kernel p{X) that 



satisfies (13 1 for all A € [0 2], or equivalently p{l) that satisfies (15 1 for all / G [—1 1]. Following Daubechies' 



approach ||T4|. we propose a maximally-flat design, in which we assign K roots to p{X) at the lowest eigenvalue 
(i.e., at A = 0). Subsequently, we select ^(A) to be the shortest length polynomial, which has K roots at A = 



and satisfies (15 1. This implies that p{l + /) has K roots at / = —1, and can be expanded as: 

k 

p{l + l) = {l + lf^r.^r. 



(24) 



m=0 



R{1) 



where R{1) is the residual k degree polynomial of /. By Corollary 3.2 there always exist such a polynomial R{1) 



On the other hand. Proposition [T] says that any p{l + /) that satisfies (15 1 has to be an odd-degree polynomial. 
Hence, p{l + /) can also be expanded as: 

M 



p(l + /) = l + J]c2„+l/2"+l. 



(25) 



n=0 



for a given M. Comparing ( |24| ) and ( [25] ), we get: 

k 



M 



(26) 



m=0 



n=0 



Comparing the constant terms in the left and right side of (26 1, we get ro = 1. Further, comparing the highest 



powers on both sides of ( 26 1 we get: 



M 



K + k-l 



(27) 
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Further, the right side in (26 1 has M constraints C2n = for n = {1,2, ...K}, and the left side in (26 1 has k 



unknowns for m = {1, 2, ...k}. In order to get a unique ■p{l + 1) that satisfies ( 15 1, we must have equal number 
of unknowns and constraints, i.e, 

M = k= ^ M = K-1. (28) 



Thus, ( |26) can be written as: 

K-1 K-1 

(1 + + r^n = 1+Y1 C2n+i/'"+\ (29) 

m=l n=0 

and K — 1 unknowns can be found uniquely, by solving a linear system of K — I equations. Note that given K, 
the length of p{l) (i.e, highest degree) is K + M = 2K — 1. As an example, we design ^(A) with K = 2 zeros at 
/ = 0. In this case p{l + /) can be written as: 

p{l + /) = (! + If {I + nl) = 1 + (n + 2)1 + (1 + 2ri)/2 + 

Since p{l + l) is an odd polynomial, the term corresponding to is zero, i.e., l + 2ri = or ri = —1/2. Therefore, 
p{l + I) is given as: 

pil + l) = + (30) 

which implies that: 

p{X) = h\3-X). (31) 

In Figure [3] we plot p{X) for various values of K, and it can be seen that by increasing K, we get a p{X) representing 
a better approximation to the ideal halfband filter. 

Note that the graph-QMF designs in ||2l were based on selecting go{X) = /io(A) hence pqmfW = ^o(-^)- Thus, 
if /io(A) is a polynomial kernel then PqmfW is the square of a polynomial, and therefore should have an even 
degree. However, as proven in Proposition [TJ any polynomial p{X) which satisfies ( 13 1 is an odd-degree polynomial. 
Therefore, ho{X) in the graph-QMF designs, cannot be an exact polynomial. 

B. Spectral factorization of half-band kernel p{X) 

Once we obtain p(A) by using the above mentioned design, we need to factor it into filter kernels /io(A) and 
go{X). Since j5(A) is a real polynomial of odd degree, it has at least one real root and all the complex roots occur 
in conjugate pairs. Since we want the two kernels to be polynomials with real coefficients, each complex conjugate 
root pair of p{X) should be assigned together to either /io(A) or go{X). While any such factorization would lead to 
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1.5 



0.5 



p{X) with K roots at A. = 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Fig. 3: The spectral distribution of p{X) with zeros at A = 



perfect reconstruction biorthogonal filterbanks, of particular interest is the design of filterbanks that are as close to 
orthogonal as possible. For this, we define a criterion based on energy preservation. In particular, we compute the 
Riesz bounds of analysis wavelet transform T^, which are the tightest lower and upper bounds, ^4 > and B < oo, 
of ||Taf||2, for any graph-signal f with ||f||2 = 1. For near-orthogonality, we require A ^ B ^ 1. The bounds 
A and B can be computed as the minimum and maximum singular values of the overall analysis side transform 
Ta of the two channel filterbank. The transform matrix consists of some rows of lowpass transform Hq (those 
corresponding to the L set) and some rows of highpass transform Hi (those corresponding to the H set) and can 
be written as: 

Ta = ^(I + J«)Ho + ^(I- J«)Hi (32) 



The singular values of are also the square roots of eigenvalues of TJTq, which can be expanded as (see lUl 
for details): 



T T 



1/2 ihl{X) + hliX))^ 



C(X) 



+ 1/2 V (/ii(A)/ii(2-A)-/io(A)/io(2-A))J^P; 



(33) 



D(A) 



where is the diagonal matrix of binary function /3. In (33 1, the term D(A) consists of product terms /io(A)/io(2— A) 



and h\{X)h\(2 — A), which are small for A away from the transition band around 1 (since these are the products 
of a low pass and a high pass kernel). Further, in the transition band when A is close to 1, the value of D{X) 
/io(A) — /ii(A) is very small compared to C(A) w /io(-^) + ^iW- Therefore, we can ignore D{X) in comparison 
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to C(A), and (33 1 can be approximately reduced to: 

t:t,«1/2 ^ ihl{X) + hliX))Fx (34) 



Thus, TjTa is a spectral transform with eigenvalues 1/2(/iq(A) + /if(A)) for A G o-{B), and the Riesz Bounds 
can be given as: 



B = ^^nv\{hl{X) + hl{\)) (35) 
We define ©, as the measure of orthogonality, given as: 

For orthogonal filterbanks © = 1. We choose filters with least dissimilar lengths, and for near orthogonal designs, 
compute © for all such possible factorizations (there are (^^^) possible choices), and choose the factorization 
with the maximum absolute value of ©. Note that the computation of ©, and hence the choice of the best solution 
depends on the exact distribution of eigenvalues A in the interval [0 2], which in turn depends on the underlying 
graph. For a graph independent design, we approximate and as the lowest and highest values, respectively, 
of 1/2(/iq(A) + h\{\)) at 100 uniformly sampled points from the continuous region [0 2], respectively, and use 
these approximations to compute ©. 

C. Unity gain compensation 

In order to avoid unnecessary growth of dynamic range in the output, it is desirable to normalize the filterbanks, 
such that the impulse responses of lowpass and highpass filters have equal (ideally unity) gain. In the graph-QMF 
case llll, the orthogonality condition ensures that all filters have unity gain. However, this is not true for the proposed 
graphBior filterbanks. Similar gain compensations have also been proposed for biorthogonal DWT filterbanks (for 
example, see ifTTIl ). The JPEG2000 standard, for example, requires the lowpass filter to have unity response for 
the DC frequency (w = 0), and the highpass filter to have unity response for the Nyquist frequency (w = tt). In 
this paper, we follow similar specifications, as given in JPEG2000, to normalize the gains of graphBior filters. In 
the bipartite graph case, uj = ^ and uj = tt correspond to the lowest (A = 0) and the highest (A = 2) magnitude 
eigenvalues, respectively. Thus, given a filter H with spectral kernel h{\), the gain factor of H is equal to |l//i(0)|, 
if H is a lowpass filter, and is equal to |l//i(2)|, if H is a highpass filter. In the filterbank implementation, a 
GC block is applied at the analysis side in each channel after filtering and downsampling, and an inverse GC 
block is applied at the synthesis side prior to filtering and upsampling. As a result, the filterbank remains perfect 
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reconstruction. 



D. Nomenclature and design of graphBior filterbanks 

The proposed biorthogonal filterbanks are specified by four parameters {k^, ki, /q, h), where ko is the number 
of roots of low pass analysis kernel /io(A) at A = 0, ki is the number of roots of low pass synthesis kernel go{X) 
at A = 0, lo is the highest degree of low pass analysis kernel /io(A), and li is the highest degree of low pass 



synthesis kernel 50(A). The other two filters, namely /ii(A) and 51(A), can be computed as in (12i. Given these 
specifications, we design p{X) = hQ{X)go{X) as a maximally fiat half band polynomial kernel with K = kQ + ki 
roots at A = 0. As a result, p{X) turns out to be a 2K — 1 degree polynomial, and we factorize it into /io(A) and 
go{X), with least dissimilar lengths (i.e., we choose Iq = K and li = K — 1). We use © as the criterion to compare 
various possible factorizations, and choose the one with the maximum value of 0. This leads to a unique design of 
biothogonal filterbanks. We term our proposed filterbanks as graphBior(ko, ki). We designed graphBior filterbanks 
for various values of (A;o, ki), and we observed that designs with ko = ki stand out, as they are close to orthogonal 
and have near- fiat pass-band responses. The lowpass and highpass analysis kernels are plotted in Figure |4j and their 
coefficients are shown in Table |Ill 



graphBior(A;o,A;i) 


filter coefficients 


ko - 6, ki - 6 


hi = [-0.3864 4.0351 -17.0630 36.5763 -39.8098 17.6477 
0] 

ho = [ 0.4352 -4.9802 23.2396 -55.4662 67.2657 -29.0402 
-13.0400 7.5253 9.5267 -4.8746 -2.0616 1.2633 1.2071] 


ko-l,ki-l 


hi = [0.3115 -3.9523 21.0540 -60.3094 98.0605 -85.9222 
31.7578 0] 

ho = [ -0.4975 6.8084 -39.6151 126.2423 -234.3683 
241.5031 -97.6557 -46.2635 62.1232 -19.3648 -2.0766 
6.5886 -4.5632 0.5775 1.5614] 


fco = 8, fci = 8 


hi = [-0.3232 4.7284 -29.7443 104.3985 -221.0705 

282.7915 -202.6283 62.8477 0000000 0] 

ho = [ 0.4470 -6.9872 47.5460 -183.6940 440.0924 

-670.0905 643.3979 -396.0713 209.9824 -154.0976 

92.8617 -30.8228 16.6112 -12.7664 3.2403 -0.0284 

1.3793] 



TABLE II: Polynomial expansion coefficients (highest degree first) of graphBior {ko, ki) filters (approximated to 
4 decimal places) on a bipartite graph. Refer to the Matlab code for more accurate coefficients. 



IV. ZERODC graphBior Filterbanks 

In many application such as images, videos and wireless sensor networks etc., where the underlying graph 
resides in a physical space, an all constant signal (a DC signal) has a physical interpretation, and the wavelet 
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graphBior(7,7) 


^hl(^) 

^hO(?.) 
^C(X) 
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graphBior(8.8) 
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-^D{X) 
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graphBior(10,10) 


^h1(>.) 

-^hO(X) 
-^C(X) 
-^D(X) 
-^p{k) + p{2-k) 
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Fig. 4: Spectral responses of graphBior(A;0) ^i) filters on a bipartite graph. In each plot, ho{X) and hi{X) are 
lowpass and highpass analysis kernels, C(A) and D{X) constitute the spectral response of the overall analysis 



filter 

Tai 'IS in (|33[). For near-orthogonality D{X) ~ and C(A) ~ 1. Finally, {p{X) +p(2 — A))/2 represents 



perfect reconstruction property as in (15 1, and should be constant equal to 1, for perfect reconstruction. 



filters are designed to be orthogonal to the dc signal. In our proposed nonzeroDC graphBior design, wavelet 
transforms (highpass transforms) are designed to be orthogonal to the eigenvector corresponding to eigenvalue 
of the normalized Laplacian matrix C, which is where D is the degree matrix. If the graph is almost 

regular (i.e., has almost same degree at all nodes), then the vector is almost constant. To obtain a zero DC 



18 



response for other non-regular graphs we propose zeroDC graphBior filterbanks designs. Construction-wise, the 
only difference between zeroDC graphBior filterbanks and nonzeroDC graphBior filterbanks is that the former are 
designed using random-walk Laplacian matrix Cr while the latter are designed using normalized Laplacian matrix 
C The two Laplacian matrices are similar, hence their eigenvalues are identical. Therefore, no change is needed 
in the design if spectral kernels designed above. In order to compute zeroDC filterbanks, we simply replace C by 
Cr in (|7]), i.e., 

Hj'j — h{(^Cr^ 

Gri = giiCr), (37) 

and replacing symmetric filters Hj and Gj by H^i and Gri, respectively, in the two-channel nonzeroDC filterbank 
implementation shown in Figure [T] The following results describe the properties of proposed zeroDC filterbanks: 
Proposition 2 (zero DC response): For any connected graph, if the spectral kernel /i(A) is such that /i(0) = 0, 
then the transform = h{Cr) has a zero DC response, i.e., H^l = 0. Further, if h{X) is a polynomial kernel 
then the transforms = h{Cr) and H = h{C) are related as: 

= D-^/2jjd1/2^ (33^ 

Proof: The random walk Laplacian matrix Cr can be diagonalized as: 

Cr = D-1/2UA(D-1/2U)-1 = D-i/2uAUTd1/2. 

Therefore, any function of Cr can be written as: 

H, = h{Cr) = D-V2u/i(A)U^Di/2 = D-i/2/i(£)Di/2 = Ti 

Further, if is an eigenvector of normalized Laplacian matrix C corresponding to eigenvalue A;, then by definition 
Hu/ = /i(A/)u;. Since uq = D^/^l is the eigenvector of C with eigenvalue 0, this implies: 

H^l = D^/^HD^/^l = Aol = 0. 



Thus Hf. has zero DC response. ■ 
Proposition 3 (perfect reconstruction property): The zeroDC filterbanks designed using graphBior spectral 
kernels are also perfect reconstruction. 



19 



Proof: Similar to ([8]l, the overall transfer function of the zeroDC filterbank; can be written as: 

f = -Gr.o(I + J/3)Hrof + -Gr.l(I — J^)Hrlf 

= -^{^roi^rO + G.rlHrl)f + -(GroJ/jHrO — GrlJ/3Hrl)f. (39) 



Using the similarity relation given in (38 1, we can simplify (39 1 as: 



f = ^(D-V2GoDV2D-i/2HoDV2^D-i/2GiDV2j3-i/2H^j3i/2)f 
+ ^(D-i/2GoDV2j^d-i/2HoD1/2 _ D-V2GiDi/2j^D-V2HiDi/2)f. 



(40) 



In (40i, the matrices D^/^, J/3, and D are diagonal matrices and hence commute with each other. Therefore, 

D^/2j^j3-i/2 ^ = 3p (41) 



Thus, (40 1, can be simplified as: 

f = ^(D-V2GoHoDi/2 + D-V2GiHiDV2)f 
+ ^(D-i/2GoJ^HoD1/2 _ D-i/2GiJ^HiDi/2)f 

= D-i/2(T,, + T,ii,,)Di/2f, (42) 

where Teg and Ta^jas correspond to the overall transfer function of nonzeroDC filterbanks, as defined in ([8]l. 
Therefore, the zeroDC filterbank implementation is equivalent to pre -multiplying the input by D~^/^ and post- 
multiplying the output by D^/^, and if the nonzeroDC filterbank is PR (i.e., T^q + T alias = cl ) then the 
corresponding zeroDC filterbank is also PR. ■ 
Proposition 4 (Riesz bounds): The zeroDC filterbanks form a Riesz basis with lower bound A^y dmin/dmax 
and upper bound B^y dmax/dmin, where A and B are the lower and upper bounds of the Riesz basis formed by 
corresponding nonzeroDC graphBior filterbanks. 

Proof: Referring to Figure [T] the wavelet coefficient vector w produced in the zeroDC filterbanks can be 
written as: 

Wr = TrJ = ^(I-J/3)H,of +^(I + J/3)H,if 

= -(H-rl + Hr.o)f + -J/3(Hrl — Hr-o)f 

= ^D-V2(Hi + Ho)DV2f + 1d-V2j^(h, - Ho)DV2f 
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This implies that the n^^ output can we written as: 

N 

dn 



Wr[n] = ^ ^I^Ta{n,m)f[m\ (44) 

m=l 



Note that if the graph is almost regular, i.e., ^ w 1, then Wr[n\ f« X]m=i "^)/[''^] ~ w[n\, where w[n] 



d 

Ah 



is the n output of the corresponding nonzeroDC filterbank. In order to obtain a worst-case bound, if we define 



fo = D^/^f, and wd = D^/^w, then (43) can be written as wd = Tafo- Thus, if the corresponding nonzeroDC 
filterbank is biorthogonal with Riesz bounds A and B, then A||fD|| < ||wd|| < i?||fD|| (the 2-norm). However, 

N N N 



dminy^^w'^ii) < ||wd|P = diw'^{i) < dmax {i) 

i=l i=l i=l 

N N N 

dmin < llfDil' = <d^„, (45) 

1=1 1=1 i=l 

where dmin is the minimum degree in the graph (1 if there is an isolated node), and dmax is the maximum degree. 



Using (45 1, we obtain: 



dmm||w|p < ||wd||^ < -B^||fD|P < -B^f^max | |f | P 
dmin^^llfll^ <^^||fD||^ < IIwdIP <d^ax||w||2, (46) 



and 

dmax J \ d' 

Thus, the zero graphBior filterbanks defines a Riesz basis in the graph-signal space, with lower bound Ar = 
Ay/ dmin /dmax and uppcr-bound Br = B^/d 

max I dmin- ^ 

Note that for regular graphs dmin = dmax, hence = {A, i?}. However, for irregular graphs the measure 

of orthogonality A,. = A^jBr = {dmin/dmax)& tend to be smaller than 0, which implies that the basis functions 
in zeroDC filterbanks are more coherent than the basis functions in nonzeroDC filterbanks. This is also confirmed 



empirically in Table III 



The decision of whether to use zeroDC graphBior or nonzeroDC graphBior filterbanks, depends upon the 
interpretation of an all constant signal 1 and its degree normalized form in the context of the problem. 

For example, in graphs arising from physical domains (sensor networks, transport networks, images and videos 
etc.), the graph signals are often nearly constant (or piecewise constant). In these cases, the all-constant signal 
should be preserved as the lowpass signal, and therefore the zeroDC filterbanks should be preferred over the 
nonzeroDC filterbanks. On the other hand, it is shown in iTSll that for highly irregular graphs (such as online social 
networks, Internet etc.) the spectral analysis gets influenced by the presence of high degree nodes, and thus misses 
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the structure around low degree nodes. Therefore, the degree normaUzed nonzeroDC filterbanks should be used for 



these cases. All the examples presented in Section VI belong to the former category (i.e., they arise in physical 



domains). Therefore, the zeroDC filterbanks are found to perform better than the nonzeroDC filterbanks. 

V. Multi-dimensional and multi-resolution implementations 

So far we have described how to implement graphBior filterbanks on bipartite graphs. This is because bipartite 
graphs provide perfect reconstruction conditions in terms of simple conditions on spectral responses in these 
filterbanks. However, not all graphs are bipartite. For arbitrary graphs, we proposed in |[T9l to decompose 
the graph G into K link-disjoint bipartite subgraphs, each defined on the entire set of vertices and their union 
covering almost all of the links in the graph. Consequently, we implemented filtering/downsampling operation in K 
stages, restricting the operations in each stage to only one bipartite graph. An example of 2-dimensional bipartite 
subgraph decomposition is shown in Figure |5j in which the graph G is divided into 4 clusters LL, LH, HL and 
HH. The first bipartite graph Bi corresponds to partitions LI = LL U LH and HI = HL U HH, and all the 
links connecting nodes in the two partitions. Subsequently, these links are removed from G and the second bipartite 
subgraph B2 corresponds to partitions L2 = HL U LL and H2 = LH U HH, and all the links between L2 and 
H2 from the remaining set of links. The remaining links are either discarded or used to further compute third and 
fourth bipartite subgraphs etc. The block diagram of a 2 "dimensional" graphBior filterbank is shown in Figure l6] 





Fig. 5: Two dimensional decomposition of a graph. 



where a dimension is interpreted as filtering and downsampling on a single bipartite subgraph. Note that this 
design is analogous to separable filterbank implementation on regular multidimensional signals. For example in the 
case of separable transforms for 2D signals, filtering in one dimension (e.g., row-wise) is followed by filtering of 
the outputs along the second dimension (column-wise). Moreover, the separable graphBior filterbanks are PR for 
any arbitrary partitions LL, LH, HL and HH induced on the graph. The choice of a specific bipartite subgraph 
decomposition depends on various factors. For highly structured graphs such as graph representation of regular 



signals, the bipartite subgraphs which preserve the structure are more useful (see, for example. Section VI-B 1. For 
arbitrary graphs, there can be various criteria. One criterion is to compute a graph decomposition that generates 
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ho 



So -\ 



Fig. 6: Separable two-dimensional filterbank on graphs. The graph is first decomposed into two bipartite subgraph 

as shown in Figure |5] The binary function /3i is such that l3i{Hi) = 1 and = —1. Similarly the binary 

function /32 is such that (3i{H2) = 1 and /3i(L2) = —1. For each bipartite graph the graph transform pair {ho, hi} 
forms the analysis low-pass and analysis high-pass graphBior filters respectively and {go,gi} axe, corresponding 
synthesis filters. GC: gain-compensation block, GC~^: inverse GC block. 



minimum number of bipartite subgraphs whose union covers all the links in the graph. An example of decomposition 

scheme based on such criterion is Horary 's algorithm proposed in ||2l, which provides a \log2K~\ bipartite subgraph 

decomposition of a K-colorable graphj^ Another criterion introduced in |[T9l . proposes subgraph decompositions 

so that the neighborhood sets of each node on different bipartite subgraphs are maximally disjoint. This leads 

to uncorrelated filtering operations on different graphs. However, whether the above mentioned decomposition 

schemes are optimal in some sense, or more generally whether there are other ways to extend graphBior filterbanks 

to arbitrary graphs, is part of our on-going research. 

The multiresolution decomposition (MR) property in graphs implies successive coarser approximations of the 

graph and graph signal. For example, in a 1 -dimensional implementation, the output samples in the set L are treated 

as signal for the next resolution level, and the vertices in L are reconnected to form a downsampled graph that 

preserves properties of the original graph such as the intrinsic geometric structure (e.g., some notion of distance 

between vertices), connectivity, graph spectral distribution, and sparsity. The graph coarsening problem has received 

a great deal of attention from graph theorists, and, in particular, from the numerical linear algebra community (see 

II20I . lIlTll and the reference therein). Further, Pesenson (e.g., ll22l ) has leveraged the analogy between the graph 

Fourier transform and the classical Fourier transform to extend the concept of bandlimited sampling to signals 

defined on graphs. Namely, certain classes of signals can be downsampled on particular subgraphs and then stably 

reconstructed from the reduced set of samples. In our designs, any of the above mentioned coarsening scheme can 

be used to compute the graph at the next level. 

'a A"-colorable graph can be divided into K clusters such that there are no links connected nodes in the same clusters. [ • ] is the ceiling 
operator. 
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VI. Experiments 

A. Performance comparisons of two-channel filterbanks on graphs 

In order to compare various graphBior designs proposed in this paper and previously proposed graphQMF designs, 
we simulate M instances of random graphs. In all the experiments the random graphs are bipartite graphs with 300 
nodes in each partition and probability of connection 2log{N) / N . The isolated vertices in the graph are removed 
in each realization. 

In order to show the trade-off between vertex domain and spectral domain localizations, we plot in Figure [7] the 
spatial spread (|2]) and spectral spread ([5]) of various two-channel spectral filterbanks on M = 10 instances of random 
bipartite graphs. We first observe that that the graph-QMF based on ideal half-band kernels (magenta diamonds in 
the plot) have very small spectral spread but very large spatial spread, as compared to other designs. This is due to 
the brick-wall spectral response of these filterbanks. The same graphQMF filterbanks when designed using smooth 
Meyer kernel based half-band filters (black squares in the plot), have lower spatial spread (though still higher 
than most of the graphBior filterbanks) but higher spectral spread. However, both of these designs do not have a 
compact support. On the other hand, the proposed graphBior filterbanks exploit the spatial/spectral tradeoff better 
and have compact support support. The filters with smaller filterlengths are spatially more localized but spectrally 
less localized than the filters with higher filterlengths. The filter length of graphBior designs is chosen to be the 
maximum of the two filter lengths (i.e, K). Among graphBior designs, the zeroDC filterbanks (red triangles in the 
plot) perform slightly worse than the nonzeroDC filterbanks, which is due to the extra normalizations introduced 
in the formers to make their DC response zero. 

The exact graphQMF filterbanks provide PR but are not compact support. In [2], we proposed polynomial 
approximation of the exact graphQMF kernels which are compact support but results in some reconstruction 
error. A comparison between proposed graphBior filterbanks and the graph-QMF filterbanks, in terms of perfect 



reconstruction error (SNR) and orthogonality (G) is shown in Table III The reconstruction SNR and orthogonality 
G are computed as an average over 20 instances of randomly generated graph-signals on M = 10 random bipartite 



graphs. It can be seen from Table III that all graphBior designs provide perfect reconstruction {SNR > lOOdB). 
The graph-QMF filters in comparison are closer to orthogonal (i.e., G almost 1), but have considerably lower 
reconstruction SNR. We now consider some applications of our proposed filterbanks. 



B. Graph based image processing 

In this section, we describe an application of proposed graphBior filterbanks for image-analysis. This is an 
extension of our previous work in ||2l, |[23l . where we proposed a graph based edge-aware representation of image- 
signals. While standard separable extensions of wavelet filterbanks to higher dimensional signals, such as 2-D 
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Fig. 7: The spatial vs. spectral spread a) of the highpass filters, and b) of the lowpass filters. The spatial/spectral 
coordinates for graphQMF filterbanks are same in both plots, since the lowpass and highpass filters are symmetric 
around A = 1. However, the lowpass and highpass filterbanks in graphBior designs are neither symmetric nor 
equal length. Therefore, the spatial/spectral spreads of the two channels are different. The dashed line is a 
quadratic polynomial fit of the data points in the least square sense. 
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4 


32.20 


0.98 


286.84 


0.88 


286.54 


0.70 


8 


32.25 


0.98 


282.89 


0.87 


282.71 


0.66 


10 


42.17 


1.00 


270.05 


0.81 


270.00 


0.65 


14 


48.09 


1.00 


230.83 


0.85 


230.73 


0.64 


16 


44.78 


0.99 


222.08 


0.94 


222.05 


0.64 


18 


45.23 


0.99 


190.53 


0.92 


190.43 


0.63 


20 


54.61 


1.00 


170.78 


0.94 


170.68 


0.63 



TABLE III: Comparison between graph-QMF filterbanks (polynomial approximations) and graphBior filterbanks 

on random bipartite graphs. 



images, provide useful multi-resolution analysis, they do not capture the intrinsic geometry of the images. For 
example, these extensions can capture only limited (mostly horizontal and vertical ) directional information. Images 
can also be viewed as graphs, by treating pixels as nodes, pixel intensities as graph-signals, and by connecting 
pixels with their neighbors in various ways. The advantage of formulating images as graphs is that different graphs 
can represent the same image, which offers flexibility of choosing the graphs that have useful properties. In 111, 
we proposed an 8-connected graph representation of images, in which each pixel is connected to 8 of its nearest 
neighbors (4 diagonal, 2 vertical and 2 horizontal) as shown in Figure [8] The graph is not bipartite, but can be 
decomposed into two bipartite subgraphs, one containing links in the horizontal and vertical direction and other in 
the diagonal directions. The proposed graphBior filterbanks can then be applied as two "dimensional" filterbanks 
as in Figure ^ The advantage of using graphBior wavelet filterbanks as against standard separable filterbanks, is 
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G B_l B_2 

Fig. 8: Two dimensional decomposition of 8-connected image-graph 



that the former provide more filtering directions (diagonal and rectangular) than the latter (only rectangular), at the 
same order of computational complexity. For multi-resolution analysis, the downsampled set of nodes in the LL 
channel, are again connected to 8 of their neighbors, to create a downsampled graph and the graphBior filterbanks 
are implemented iteratively on the downsampled graphs. Thus the downsampling ratio at each level is same in both 
graphBior filterbanks and standard separable filterbanks. 

In |[23l . we proposed an edge-aware implementation for piece-wise smooth images, in which the bipartite 
subgraphs obtained in Figure [8} can be simplified by removing the links between pixels across which the pixel 
intensity changes drastically. These links can be found using any standard edge-detection algorithm (we use Canny 
edge detection in our experiments and remove connected components less than 50 pixels before computing the 
graph). The advantage of edge-aware graph representations is that it avoids filtering across edges, which leads to a 
very significant reduction in the number of large coefficients near edge (and thus corresponding reductions in rate). 
Note that in a compression application, this would require generating an edge map at the encoder and then sending 
it to the decoder. However, recent work Il24l . ||25l using transforms based on similar edge-map information have 
shown that even with the extra overhead of sending the edge map we can achieve reductions in overall transmitted 
rate. 

In order to demonstrate the advantage of graph based implementation of proposed filterbanks, we choose coins.png 



image as shown in Figure 9(a) with many round shaped coins. We implement graphBior filterbanks of length 10, 
(i.e., graphBior{5, 5)), and compare them against standard separable CDF 9/7 wavelet filterbanks in a non-linear 
approximation of images, using 4 resolution levels. Figure |9] shows the reconstruction of coins.png using all lowpass 
coefficients and top 4% of the highpass coefficients in terms of magnitude. Since the standard separable wavelet 
filterbanks filter only in horizontal and vertical directions, they produce lots of large magnitude wavelet coefficients 
(and hence blurring artifacts) near the edges (see Figure 9(b)[ ). The zeroDC graphBior filterbank implementation 



on the regular 8-connected image graph (Figure 9(c) I does slightly better since it also provide filtering in diagonal 
directions. However, the best performance in terms of reconstruction quality is observed for proposed edge-aware 



zeroDC graphBior filterbanks, especially in preserving the edge structure (see Figure 9(d) i. This is due to the fact 
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(d) (e) (f) 

Fig. 9: Reconstruction of "Coin.png" (512 x 512) from all lowpass coefficients and 3% highpass coefficients after 

a 4-level decomposition, (a) Original image, (b) standard CDF 9/7 filters, (c) zeroDC filterbanks on regular 
8-connected image graph (d) zeroDC filterbanks on edge-aware image graph (e) nonzeroDC filterbanks on regular 
8 connected image graph and (f) nonzeroDC filterbanks on edge-aware image graph. 



that the underlying graphs in this approach are disconnected at the edges, and hence the filtering operations do not 
cross the edges. Theoretically, the nonzeroDC filterbanks should perform almost the same as zeroDC filterbanks 
for regular degree graphs. The 8-connected image-graphs are almost regular except at the boundaries, and edges, 
and we observe in Figures 9(e) and 9(f)[ that significant ringing artifacts are produced near these places, when 
using nonzeroDC filterbanks. The problem of boundary artifacts also arises when using standard filterbanks on 
images, which is usally solved by providing signal extensions at the boundaries. Whether such signal extensions 



can be proposed for graph representation of images, is an open issue. Figure 10 shows PSNR and SSIM 11261 
values plotted against fraction of detail coefficients used in the reconstruction of coins. png image, and it can be 
seen from both the plots that zeroDC graphBior filterbanks perform better (up to 2dB better in PSNR) than the 
standard CDF 9/7 filterbanks. Thus, the results show that the proposed graphBior filterbanks provide advantages 
over the standard wavelet transforms, with the same order of computational complexity. 



C. Compression and Learning on arbitrary graphs 

The proposed filterbanks are useful in analyzing and compressing signals defined on arbitrary graphs. As a proof 
of concept, we implement proposed graphBior filterbanks on the Minnesota traffic graph used in [ZJ. The graph 
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(a) (b) 

Fig. 10: Reconstruction of coins.png image from all low-pass coefficient and a fraction of wavelet coefficients 
(sorted in the order of magnitudes). The fraction value is plotted on the x-axis. (a) PSNR of the reconstructed 

images, (b) SSIM of the reconstructed image. 



is shown in Figure [TTJa), and the graph signal to be analyzed is shown in Figure [TTJb), where the color of a 
node represents the signal value at that node. The graph is perfectly 3-colorable and hence, it can be decomposed 



using Harary's decomposition Q into [^052 (3)] = 2 bipartite subgraphs, which are shown in Figure 11 c-d), and 
a 2-dimensional graphBior filterbank given in Figure [6] with filterlength = 10 is implemented on the graph. 

The output in the 4 channels can be interpreted as follows: the LL channel providing a smooth approximation 
of the original signal on a subset of nodes, and the remaining channels providing details required for perfect 
reconstruction. Moreover, the total number of outputs in all channels is equal to the total number of input samples, 
hence the transform is critically sampled. The HL channel does not sample any output and is empty. The output 



coefficients of both zeroDC and nonzeroDC filterbanks in LL, LH and HH channels are shown in Figure 12 



Note that the graph-signal is piece-wise constant, hence the proposed zeroDC filterbanks (bottom row in Figure 12 1 



provide a sparser approximation than the nonzeroDC filterbanks (top row in Figure 12 1. As a result, the non-linear 



approximation of the graph signal with only 1% highpass coefficients (and all low pass coefficients) provide better 



SNR when using zeroDC filterbanks than when using nonzeroDC filterbanks as shown in Figure 13 



VII. Conclusions 

In this paper we have presented novel graph-wavelet filterbanks that provide a critically sampled representation 
with compactly supported basis functions. The filterbanks come in two flavors: a) nonzeroDC filterbanks, and 
b) zeroDC filterbanks. The former filterbanks are designed as polynomials of the normalized graph Laplacian 
matrix, and the latter filterbanks are extensions of the former to provide a zero response by the highpass operators. 
Preliminary results showed that the filterbanks are useful not only for arbitrary graph but also to the standard 
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Input Signal 




(c) (d) 



Fig. 11: (a) The Minnesota traffic graph G, and (b) the graph-signal to be analyzed. The colors of the nodes 
represent the sample values. (c)(d) bipartite decomposition of G into two bipartite subgraphs using Harary's 

decomposition. 

regular signal processing domains. Extensions of this work will focus on the application of these filters to different 
scenarios, including, for example, social network analysis, sensor networks etc. 
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